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Following the trichromatic sub-Doppler cooling scheme on the D\ line of 40 K recently demon- 
strated in D. R. Fernandes et al. (2012), we introduce a similar technique for 7 Li atoms and obtain 
temperatures of 60 /iK while capturing all of the 5 x 10 8 atoms present from the previous stage. 
We investigate the influence of the detuning between the the two cooling frequencies and observe 
a threefold decrease of the temperature when the Raman condition is fulfllled. We interpret this 
effect as arising from extra cooling due to long-lived coherences between hyperfine states. Solving 
the optical Bloch equations for a simplified, A-type three-level system we identify the presence of an 
efficient cooling force near the Raman condition. After transfer into a quadrupole magnetic trap, 
we measure a phase space density of ~ 10~ 5 . This laser cooling offers a promising route for fast 
evaporation of lithium atoms to quantum degeneracy in optical or magnetic traps. 

PACS numbers: 37.10.De, 32.80. Wr, 67.85.-d 
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Introduction 

Lithium is enjoying widespread popularity in the cold- 
atom trapping community thanks to the tunability of 
its two-body interactions and its lightness. Both the 
fermionic and bosonic isotopes of lithium feature broad, 
magnetically-tunable Feshbach resonances in a number 
of hyperfine states [1 J. The presence of these broad reso- 
nances makes lithium an attractive candidate for studies 
of both the Fermi- and Bose-Hubbard models [2 and 
the strongly correlated regime for bulk dilute gases of 
Fermi [3 or Bose [3H6] character. Its small mass and 
correspondingly large photon-recoil energy are favorable 
factors for large area atom interferometers [7] and pre- 
cision frequency measurements of the recoil energy and 
fine structure constant [8 j. Under the tight-binding lat- 
tice model, lithium's large photon-recoil energy leads to 
a larger tunneling rate and faster timescale for super- 
exchange processes, allowing for easier access to spin- 
dominated regimes [9 J. Finally, lithium's small mass re- 
duces the heating due to non-adiabatic parts of the col- 
lision between ultracold atoms and Paul-trapped ions. 
This feature, together with Pauli suppression of atom- 
ion three-body recombination events involving 6 Li [TO] . 
potentially allows one to reach the s- wave regime of ion- 
atom collisions [TT] . 

However, lithium, like potassium, is harder to cool us- 
ing optical transitions than the other alkali atoms. The 
excited-state structure of the D^ transition in lithium 
lacks the separation between hyperfine states for stan- 
dard sub-Doppler cooling techniques such as polarization 
gradient cooling [T2HI4] to work efficiently. Recently, 
it has been shown by the Rice group that cooling on 
the narrow 2Si/2 — >• 3P3/2 transition produces lithium 
clouds near 60 /iK, about half the Z^-hne Doppler cooling 
limit [15] , and can be used for fast all-optical production 
of a 6 Li quantum degenerate Fermi gas. However, this 
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FIG. 1: (Color online) The Di-line for 7 Li. The cooling 
scheme has a strong coupling laser (principal beam, black 
solid arrow) 62 blue detuned from the \F = 2) — » \F' — 2) 
transition and a weak coupling laser (repumper, gray solid 
arrow) 5\ blue detuned from the \F = 1) — >> \F' — 2) transi- 
tion. The repumper is generated from the principal beam 
by an electro-optical modulator operating at a frequency 
803.5 + 6/2tv MHz where 6 = 61 - 6 2 . 



approach requires special optics and a coherent source 
at 323 nm, a wavelength range where power is still lim- 
ited. Another route is to use the three-level structure 
of the atom as implemented previously in neutral atoms 
and trapped ions [T6H22] . The three- level structure of- 
fers the possibility of using dark states to achieve tem- 
peratures below the standard Doppler limit, as evidenced 
by the use of velocity-selective coherent population trap- 
ping (VSCPT) to produce atomic clouds with sub-recoil 
temperatures [23 j. In another application, electromag- 
netically induced transparency has been used to demon- 
strate robust cooling of a single ion to its motional ground 
state [T91I24]. 

In this paper, we implement three-dimensional bichro- 



matic sub-Doppler laser cooling of 7 Li atoms on the 
D\ transition. Figjl] presents the 7 Li level scheme and 
the detunings of the two cooling lasers that are ap- 
plied to the atoms after the magneto-optical trapping 
phase. Our method combines a gray molasses cooling 
scheme on the \F — 2) — > \F' = 2) transition citewei- 
demullerl994,fernandes2012 with phase-coherent ad- 
dressing of the \F — 1) — >• \F' = 2) transition, creating 
VSCPT-like dark states at the two-photon resonance. In- 
stead of UV laser sources, the method uses laser light that 
is conveniently produced at 671 nm by semi-conductor 
laser sources or solid-state lasers [25j [26] with sufficient 
power. This enables us to capture all of the ~ 5 x 10 8 
atoms from a MOT and cool them to 60 fiK in a duration 
of 2 ms. 

We investigate the influence of the relative detuning 
between the two cooling lasers and observe a threefold 
decrease of the temperature in a narrow frequency range 
around the exact Raman condition. We show that extra 
cooling arises due to long-lived coherences between hy- 
perfine states. We develop a simple theoretical model 
for a sub-Doppler cooling mechanism which occurs in 
atoms with a A-type three-level structure, in this case, 
the F = 1, F = 2, and F' = 2 manifolds of the D 1 
transition in 7 Li. The main physical cooling mechanism 
is contained in a ID bichromatic lattice model. We first 
give a perturbative solution to the model and then verify 
the validity of this approach with a continued fraction 
solution to the optical Bloch equations (OBEs). 



I. EXPERIMENT 

The stage preceding D\ sub-Doppler cooling is a com- 
pressed magneto-optical trap (CMOT) in which, starting 
from a standard MOT optimized for total atom number, 
the frequency of the cooling laser is quickly brought close 
to resonance while the repumping laser intensity is di- 
minished in order to increase the sample's phase space 
density [27 . The CMOT delivers 5 x 10 8 7 Li atoms at 
a temperature of 600 jaK. The atoms are distributed 
throughout the F — 1 manifold in a spatial volume of 
800- /im 1/e width. Before starting our D 1 molasses cool- 
ing, we wait 200 /is to allow any transient magnetic fields 
to decay to below 0.1 G. The light used for Di cool- 
ing is generated by a solid-state laser presented in [25] . 
The laser is locked at frequency c^2, detuned from the 
\F = 2) -)► \F' = 2) D 1 transition in 7 Li by 5 2 . It is then 
sent through a resonant electro optical modulator (EOM) 
operating at a frequency near the hyperfine splitting in 
7 Li, z^eom = 803.5 MHz -\-S/2tt. This generates a small 
amplitude sideband, typically a few % of the carrier, at 
frequency uji . We define the detuning of this frequency 
from the \F = 1) — >• \F' = 2) transition as S\ (such that 
S = Ji — #2), as shown in Fig. [I] Using about 150 mW of 
671 nm light we perform a three-dimensional D\ molasses 
as in [28 , with 3 pairs of <r + — o~ counter-propagating 
beams. The beams are of 3.4-mm waist and the intensity 
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FIG. 2: (Color online) a)Typical temperature of the cloud as a 
function of the repumper detuning for a fixed principal beam 
detuned at £1 = 4.5r = 2tt x 26.4 MHz. The dashed verti- 
cal line indicates the position of the resonance with transition 
\F = 2) — >- \F' = 2), the dotted horizontal line shows the typi- 
cal temperature of a MOT. b) Zoom-in on the region near the 
Raman condition with well-aligned cooling beams and zeroed 
magnetic offset fields, c) Minimum cloud temperature as a 
function of repumper power. 



(/) of each beam is / > 45/ sat where I sat = 2.54mW/cm 2 
is the saturation intensity of the D 2 cycling transition in 
Lithium. We capture all of the atoms present after the 
CMOT stage into the D\ gray molasses. The 1/e lifetime 
of atoms in the molasses is >50 ms. After being cooled 
for 1.5-2.0 ms, the temperature is a low as 40 /uK without 
optical pumping or 60 jaK after optical pumping into the 
\F = 2 1 rriF — 2) state for imaging and subsequent mag- 
netic trapping. In contrast with [28], we find no further 
reduction in the steady-state temperature by slowly low- 
ering the light intensities after the initial 2.0 ms. 

During the molasses phase, we find a very weak depen- 
dence on the principal laser detuning for 3r < S 2 < 6I\ 
For the remainder of this article, we will use a prin- 
cipal laser detuning of S 2 = 4.51" = 2i\ x 26.4 MHz. 
In Fig. [2|a), the temperature dependence upon the re- 
pumper detuning is displayed for typical conditions. For 
— 9 < S/T < —6, the temperature drops from 600 fiK (the 
CMOT temperature) to 200 fiK as gray molasses cooling 
gains in efficiency when the weak repumper comes closer 
to resonance. For —6 < S/T < — 1, the cloud temperature 
stays essentially constant but, in a narrow range near the 
position of the exact Raman condition (5 = 0), one no- 
tices a sharp drop of the temperature. For S slightly blue 
of the Raman condition, a strong heating of the cloud oc- 
curs, accompanied by a sharp decrease in the number of 
cooled atoms. Finally for 6 > T, the temperature drops 
again to a level much below the initial MOT temperature 
until the repumper detuning becomes too large to pro- 
duce significant cooling below the CMOT temperature. 

Figs, ^[b) and (c) show the sensitivity of the temper- 



ature minimum to repumper deviation from the Raman 
condition and repumper power, respectively. The tem- 
perature reaches 60 jiK in a ±500 kHz interval around 
the Raman resonance condition. After taking the data 
for Fig. [2|a), the magnetic field zeroing and beam align- 
ment were improved, which accounts for the frequency 
offset and higher temperature shown in (a) relative to 
(b) and (c). The strong influence of the repumper around 
the Raman condition with a sudden change from cooling 
to heating for small and positive Raman detunings moti- 
vated the study of the bichromatic-lattice effects induced 
by the A-type level configuration which is presented in 
the next section. 



II. MODEL FOR HYPERFINE RAMAN 

COHERENCE EFFECTS ON THE COOLING 

EFFICIENCY 

In order to understand how the addition of the sec- 
ond manifold of ground states modifies the gray molasses 
scheme, we analyze a one-dimensional model based on 
a A-type three-level system schematically represented in 
Fig§ 




FIG. 3: The A level scheme. An intense standing wave with 
Rabi frequency Q2, and a weaker standing wave with Rabi fre- 
quency Qi, detuning 5i, illuminate an atom with three levels 
in a A configuration. 



1 3) with a Rabi frequency ^ 2 = Y^/l/2I sat where / is 
the laser light intensity and I sat the saturation intensity 
on this transition. The repumper transition is labeled 1, 
between states |1) and |3) with Rabi frequency Q\ much 
smaller than ^2- 

The corresponding Hamiltonian for the light-atom in- 
teraction in the rotating wave approximation (at ui) is 



H a .i. = htt 2 cos(kz) (|2) (3| + h.c.) 
+ ftfticos(fc2 + 0)(|l)(3| H 
+ M 2 |2)(2| + Mi|l)(l|. 



h.c.) 



(1) 



The usual formalism used to compute the atom's dy- 
namics is to consider the light force as a Langevin force. 
Its mean value is ^(v), and the fluctuations around this 
mean will give rise to diffusion in momentum space, char- 
acterized by the diffusion coefficient V p (v) > 0. In order 
to calculate an equilibrium temperature, one needs J-(v) 
and V p (v). In the limit of small velocities the force reads 



Hv) 



■av, 



(2) 



with a the friction coefficient. If a > the force is a cool- 
ing force, in the opposite case it produces heating. For 
a cooling force the limiting temperature in this regime is 
given by 



k B T ~ V p (0)/a. 



(3) 



However, since our model is a gross simplification of 
the physical system, we do not expect to be able to quan- 
titatively predict a steady-state temperature. Instead, in 
order to reveal the physical mechanisms in action, we 
only calculate the force F{v) and the excited state popu- 
lation P33. Restricting our analysis to the force and pho- 
ton scattering rate, Y • P33, suffices to determine whether 
the action of the weak repumper serves to heat or cool 
the atomic ensemble. 

From dTl) the mean light force on the atoms is computed 
by taking the quantum average of the gradient of the 



potential, F 
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with p the 



density matrix, yielding the wavelength-averaged force T 



T{v) 



k 

2^ 



dz F(z, v) 



(4) 



A. The model 

This model includes only the F = 1, 2 hyperfine ground 
states and the F' = 2 excited state ignoring the Zeeman 
degeneracy, hence standard gray molasses cooling [29] 
will not appear in this model. The states are addressed 
by two standing waves with nearly the same frequency 
uoi ~ UJ2 — uj = kc but spatially shifted by a phase <\>. The 
principal cooling transition F = 2 — >• F' = 2 is labeled 
here and below as transition 2, between states |2) and 
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dzs'm(kz) (^2Rep23 + ^iRepi2) (5) 



The spontaneous emission rate averaged over the stand- 
ing wave is simply given by the linewidth of the excited 
state multiplied by its population: 



r = -y o d£ r to 



(6) 



So, both the force and the spontaneous emission rate are 
functions of the density matrix p, the evolution of which 



is given by the optical Bloch equations (OBEs) 
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i—p 

dr 



Hal,P 



dp\ 



l[ ti) (7) 

/ spont. emis. 



As we will be focusing on the sub-Doppler regime, we 
assume 



v<^T/k 



(8) 



with v being the velocity. The inequality holds for T <C 
13 mK for lithium. This inequality allows us to replace 
the full time derivative in the l.h.s. of (|7|) by a partial 
spatial derivative times the atomic velocity 
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dt 



d_ 
dz 



Using the notation Qi(z) = QiCos(z + ^), and setting 
h = k = 1 from here on: 



0p 2 2 

1 dz 
dpn 
1 dz 

dp 2 3 

1 dz 
dpi3 
1 dz 



dz 



-2iQ 2 (z) Im(p 2 3) + i^P33 

r 

-2iQi(z) Im(pi 3 ) + i~P33 



(9) 
(10) 



($2 - iy)p23 + H 2 (^) (p33 - P22) - Ql{z)p2l (11) 

r 

(£l - ^-)P13 + ^l(z) (p33 ~ Pll) ~ ^2(z)p 12 (12) 



(S 2 - Si)p21 + ^2(z)p 3 l - 9, 2 {z)p2 



(13) 



The solution of these equations yields the expression of 
J~(v) and r' . This semi-classical model is valid only for 
velocities above the recoil velocity v Tec = hk/m (corre- 
sponding to a temperature mv rec /kB of about 6 pK for 
Lithium). Different theoretical studies [17, 18, 20, 22, 
l30l l3T] as well as experiments [161 [32] have been per- 
formed on such a A-configuration in standing waves or 
similar systems. However in our 7 Li experiment, we have 
the situation in which the A configuration is coupled to 
a gray molasses scheme which involves a different set of 
dark states. This fixes the laser light parameters to val- 
ues that have not been studied before and motivates our 
theoretical exploration. Thus we will concentrate on the 
situation corresponding to the conditions of our experi- 
ment. 

To solve the OBEs (9|13), we first introduce a pertur- 
bative approach that enables us to point out the relevant 
physical mechanisms. We further extend the analysis by 
an exact approach in terms of continued fractions. 



B. Perturbative approach 

In our perturbative approach we choose a Rabi fre- 
quency ^2 between 2r and 41" and Qi <C T, ^2,^2 as 
the ratio of the repumper to principal laser power is very 
small, typically (f^i/f^) 2 < .03, under our experimental 



conditions. We further simplify the approach by consid- 
ering only the in-phase situation <j> = 0; any finite phase 
would lead to divergencies of the perturbative approach 
at the nodes of wave 1. The validity of these assumptions 
will be discussed in section HTCl 
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FIG. 4: (Color online) Comparison of experimental data with 
the perturbative approach results for a detuning of the pump 
82 = 27T x 26.4 MHz = 4.5r, (a) temperature versus repumper 
detuning, experiment, we indicate the MOT temparature by 
the dotted line, (b) and (c) show respectively the friction co- 
efficient a and photon scattering rate T' for 0,2 = 3.4I 1 , red 
dashed curve and 2. IT, blue solid curve. The intensity ra- 
tio (Q1/Q2) 2 is 0.02. The vertical dashed line indicates the 
position of Si =0. 



We perform an expansion in powers of the Rabi fre- 
quency ^1 and the atomic velocity, such that the com- 
plete expansion reads: 



^ = £/€W>)'- 



(14) 



n.l 



This expansion of p allows us to recursively solve the 
OBEs. Using an expansion similar to eq. (14) for the 
force, we find: 



a = -J2 jr{nA) (tt 1 ) n . 



(15) 



n=0 



We plug the perturbative solution of the OBEs into 
eq. ([5| and find, to the lowest order (n = 2) in Q±: 

- - lf|) ■ / " dzsin(z) (sizRep^ + itep& 1) )(16) 



2tt 



The spontaneous emission rate to lowest order in v and 
Qi reads 



r' 



2tt 



Jo 



j (2,0) 



(17) 



Figure [4] presents the results from (15, 17) compared 
with the experimental data. It shows that indeed a nar- 
row cooling force appears near the Raman resonance con- 
dition and that the photon scattering rate vanishes at 
exact resonance, hinting towards an increase of cooling 
efficiency with respect to the gray molasses Sisyphus cool- 
ing mechanism which achieves a temperature near 200 fiK 
over a broad range. The strong heating peak for small, 
positive repumper detuning is also a consequence of the 
negative value of a, and the heating peak shifts towards 
higher frequency and broadens for larger intensities of 
the principal laser. In contrast, the friction coefficient 
and scattering rate in the range —6 < S/T < — 3, which 
correspond to a repumper near resonance, do not seem 
to significantly affect the measured temperature. 

To gain further physical insight into this cooling near 
the Raman condition, it is useful to work in the dressed- 
atom picture. Given the weak repumping intensity, we 
first ignore its effect and consider only the dressing of 
the states |2) and |3) by the strong pump with Rabi fre- 
quency 0,2- This dressing gives rise to an Autler-Townes 
doublet structure which follows the spatial modulation 
of the standing wave: 



|2'> 

13') 



oc 



oc 



\2)-i£l 2 (z)/6 2 \3) 
-iQ 2 (z)/S 2 |2> + |3> 



(18) 
(19) 



Since the pump is relatively far detuned (in the con- 
ditions of FigJ4j^ 2 /^2 < 0.45), the broad state |3') car- 
ries little 1 2) cnaracter. Conversely the narrow state |2') 
is mostly state |2). It follows that |3') has a lifetime 
pi 3 '} ~ p while 1 2') is relatively long-lived with a spa- 
tially dependent linewidth rl 2 ^ = T {^{z) / 52) 2 which 
is always < T/6 for the parameters chosen here. In 
order to reintroduce the effects of the repumping ra- 
diation, we note that the position in 5 of the broad 
state is S^ 3 ^ c± — 62 — ^2(^) 2 /^2 and the narrow state 
S^ 2 '} ~ ^2(^) 2 /^2- As coherent population transfer be- 
tween |1) and 1 2') does not change the ensemble temper- 
ature, we consider only events which couple atoms out of 
1 2') to |1) through spontaneous decay and therefore scale 
with T| 2 /) . The rates of coupling from |1) into the dressed 
states can be approximated by the two-level absorption 
rates: 



7|i>-H2'> 



7|i>-H3') 



^{zf 



Y\ 2 '\z) 



{Y\ 2 '){z)/2Y + (5-5\ 2 ')(z)) 2 



ni(z) 



{T/2y + (5-5\*'){z)y 



(20) 



(21) 



Finally, these results are valid only in the limit \5\ > 
Tti^/fi^ ( see e -9- [33J when state |1) is weakly coupled 
to the radiative cascade. Near the Raman resonance, the 
dressed state family contains a dark state which bears 



an infinite lifetime under the assumptions made in this 
section but is in reality limited by off-resonant excitations 
and motional coupling. This dark state reads: 



|NC) = (fi 2 |l>-fii|2))/V«? + fii, 
which we must add in by hand. 



H, n> 




H, n) 




(22) 
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FIG. 5: (Color online) The cascade of levels dressed by tran- 
sition 2 with a schematical representation of state |1). Traces 
show typical cycles of atoms pumped from |1) and back de- 
pending on the detuning of wave 1. The detuning of the 
repumper modulates the entry point into the cascade of the 
dressed states, leading either a) heating or b) cooling pro- 
cesses. 

Using this toy model, we now explain the features of 
Fig. [4] and Fig. |2j Figure [5] represents the cascade of 
dressed levels where each doublet is separated by one 
pump photon. It gives rise, for example, to the well- 
known Mollow triplet. Condition (M) states that if an 
atom falls in state |3') it will rapidly decay to |2') without 
traveling a significant distance. However, the atom will 
remain in \2') long enough to sample the spatial variation 
of the standing wave and gain or lose energy depending 



on the difference of light shift between the entry and de- 
parture points, as in most sub-Doppler cooling schemes. 

Let us first analyze the spontaneous emission rate 
shown in Fig. W[c) . It reaches two maxima, the first one 
for S ~ Jl 3 ^ and the second one for S ~ S^\ and it goes 
to exactly zero at 5 = 0. The two maxima are simply 
due to scattering off the states |2') and |3'). At 5 = 0, 
I" goes to zero due to coherent population trapping in 
|NC). It is the presence of this dark state which leads to 
the reduced scattering rate of photons around 6 = and 
the suppression of the final temperature of the gas in the 
region around the Raman condition. 

The friction coefficient, Fig.[4Fb), displays a more com- 
plicated structure with variations in S. It shows a dis- 
persive shape around J' 3 ^; remains positive in the range 
Jl 3 '^ < S < 0; diverges at S = 0; and reaches negative 
values for 6 > up to 6\ 2 '\ where it drops to negligi- 
ble values. This structure for a can be explained using 
our toy model. Let us consider the different scenarios 
corresponding to both sides o f 6 near 0, they follow for- 
mally from equations ( 2Q|21 ) and the spatially varying 
linewidthof |2'). 

For the case of the repumper tuned slightly blue of the 
narrow doublet state, S > S^\ shown in Fig. Kna), the 
atoms are pumped directly from |1) into |2'). However, 
this pumping happens preferentially at the antinodes of 
the standing wave as the repumper intensity is greatest, 
the linewidth of |2 7 ) is the largest, and the light shift min- 
imizes the detuning of the repumper from the |1) —> \2 f ) 
transition for the = case considered here. On aver- 
age, the atoms exit this state at a point with a smaller 
light shift through a spontaneous emission process either 
into the cascade of dressed states or directly back to |1). 
As a result, we expect heating and a < in this region. 

For repumper detunings between 5^'^ and 0, Fig.J5[b), 
we predict cooling. For this region, the atoms are ini- 
tially pumped into |3'). Here the light shift modifies 
the relative detuning, favoring coupling near the nodes 
of the light. Spontaneous decay drops the atoms near 
to the nodes of the longer-lived |2'), and they travel up 
the potential hill into regions of larger light shift before 
decaying, yielding cooling and a positive a. These sign 
changes of a and the decreased scattering rate due to 
|NC) in the vicinity of the Raman condition explain the 
features of our perturbative model. 

We conclude this section by stating that the exper- 
imentally observed change of sign of the force close to 
the Raman condition is well described in our perturba- 
tive model. The model further reveals the importance 
of Raman coherence and the existence of a dark state. 
The dark state together with the friction coefficient as- 
sociated with cycles represented in trace 5(b) correspond 
to a cooling mechanism analogous to that of gray mo- 
lasses. In this way, the bichromatic system provides an 
additional gray molasses scheme involving both hyper- 
fine states which complements the gray molasses cooling 
scheme on the principal transition. On the other hand, 
when the friction coefficient is negative in the vicinity of 



the two-photon resonance, it turns into a heating mech- 
anism that overcomes the standard gray molasses oper- 
ating on the F = 2 —> F' = 2 transition. 

The perturbative approach successfully revealed the 
mechanisms giving rise to the experimentally observed 
additional cooling. However, it also possesses some short- 
comings. First, the divergence of a at 5 = is not 
physical; the assumption that Vt\ is the smallest scale 
in the problem breaks down when <S — >• 0. Alternatively, 
it can be seen as the failure of our model based on non- 
degenerate perturbative theory in the region where |1) 
and 1 2) become degenerate when dressed with oj\ and co>2, 
respectively. Secondly, we have only addressed the = 
case. Since the experiment was done in 3 dimensions with 
3 pairs of counter-propagating beams, the relative phase 
between the two frequencies varies spatially, and we must 
test if the picture derived at = holds when averaging 
over all phases. In order to address these limitations and 
confirm the predictions of the perturbative approach, we 
now present a continued fractions solution to the OBEs 
which does not rely on Qi being a small parameter. 



C. Continued fractions approach 

The limitations listed above can be addressed by using 
a more general approach, namely an expansion of the 
density matrix in Fourier harmonics: 



Pij 



z=+oo 
V^ (n) „ 






-% 3 



(23) 



-oo 



Injecting this expansion in (9][l3) yields recursive rela- 
tions between different Fourier components of p. Koza- 
chiov et al. [I7J[30] express the solutions of these relations 
for a generalized A system in terms of continued fractions. 
Here we use their results to numerically solve the Bloch 
equations. We then compute the force F(v) to arbitrary 
order of Qi and extract a by means of a linear fit to the 
small-v region. We then compute F(v ) and the photon 
scattering rate V averaged over the phase between the 
two standing waves. 

Fig. [6] compares a (S) obtained through the continued 
fractions approach with the results of the perturbative 
expansion for the = case. The continued fractions 
approach has removed the divergence at S = and a 
crosses zero linearly. The overall friction coefficient is 
reduced but the two methods show qualitative agreement 
in the range of S considered. At the Raman condition 
the interaction with light is canceled due to the presence 
of |NC), thus the diffusion coefficient V p in momentum 
space also cancels. To lowest order, the diffusion and 
friction coefficients scale as 



■D p * S 2 
a ~ 5 

according to (|3| the temperature scales as: 

T~5. 



(24) 
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FIG. 6: Comparison of results using the perturbative calcu- 
lation (dashed), and the continued fractions (solid) for the 
= case, with the same parameters as in Fig. Q and 

0,2 = 2.ir. 



Through this qualitative scaling argument, we show that 
even though the light action on the atoms is suppressed 
when approaching the Raman condition, we expect that 
the temperature will drop when approaching from the 
S < side, completing the physical picture derived in 
the previous section. 

Next, we analyze how a randomized phase between 
the repumping and principal standing waves, </>, modifies 
J~(y). In order to take this into account, we calculate the 
phase- averaged force: 



(Hv))^ 



i r 2n 



T{v, 4>) d(j>. 



(27) 



In Fig. [7| the phase-averaged force is plotted for 
various detunings near the Raman condition. It can be 
seen that a cooling force is present for small detunings, 
qualitatively in agreement with our perturbative model 
and with the experimental data. The force, however, 
changes sign to heating for small blue detuning, close 
to S — 0.6 T, also in qualitative agreement with the 
experimental data. We note that the cooling slope very 
close to zero velocity in the S = 0.8 T plot corresponds 
to a velocity on the order of or below the single-photon 
recoil velocity, i.e. is non-physical. 

Finally, for the </> 7^ case, |NC) varies in space and the 
motion of the atoms can couple atoms out of |NC) even 
at the Raman condition. In Fig. [8] we verify that the rate 
of photon scattering retains a minimum near the 6 = 
region after averaging over <j) by plotting (T f ) , = T (^33)0 
calculated with the continued fractions approach. Over- 
all, the friction coefficient a and photon scattering rate V 
confirm the existence of a cooling force associated with a 
decrease in photon scattering in the vicinity of the Raman 
condition for the 1-D bichromatic standing wave model. 
Thus, the continued fractions calculation has confirmed 
the physical mechanisms revealed by the perturbative ex- 
pansion and that the lowest temperatures should be ex- 
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FIG. 7: (J 7 ), in units of 1/hkT as a function of v for different 
values of 5 around 5 — 0. The horizontal scale is in units of 
the thermal velocity at T = 200 fiK, 
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FIG. 8: Continued fractions solution of the photon scattering 
rate T' = V P33 averaged over all relative phases of the re- 
pumper and principal standing waves as a function of the two- 
photon detuning 6. Velocity-dependent effects are taken into 
account here by computing an average of (r ; ) , (v) weighed 
by a Maxwell- Boltzmann velocity distribution at 200 /xK. 



pected close to S = 0, as seen in the experiment. 



III. CONCLUSION 

In this study, using bichromatic laser light near 670- 
nm, we have demonstrated sub-Doppler cooling of 7 Li 
atoms down to 60 /iK with near unity capture efficiency 
from a magneto-optical trap. Solving the optical Bloch 
equations for a simplified A level structure, we have an- 



alyzed the detuning dependence of the cooling force and 
photon scattering rate. Our analysis shows that the low- 
est temperatures are expected for a detuning of the re- 
pumping light near the Raman condition, in agreement 
with our measurements. There the A-configuration adds 
a new set of long-lived dark states that strongly en- 
hance the cooling efficiency. For T Li, this addition re- 
sults in a threefold reduction of the steady-state tem- 
perature in comparison with an incoherently repumped 
gray molasses scheme. This atomic cloud at 60/iK is an 
ideal starting point for direct loading into a dipole trap, 
where one of the broad Feshbach resonances in the lowest- 
energy states of 7 Li or 6 Li could be used to efficiently cool 
the atoms to quantum degeneracy [TSJ |34] . Alternatively, 
when the atoms are loaded into a quadrupole magnetic 
trap, we measure a phase space density of ~ 10 -5 . This 
A-enhanced sub-Doppler-cooling in a D\ gray molasses is 



general and should occur in all alkalis. Notably, we have 
observed its signature in a number of the alkali isotopes 
not amenable to polarization gradient cooling: 7 Li (this 
work), 40 K [28], and 6 Li [35]. 
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